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ABSTRACT 
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. The zodiacal light is the dominant source of the mid-infrared sky brightness seen 

, from Earth, and exozodiacal light is the dominant emission from planetary and debris 

Q I systems around other stars. We observed the zodiacal light spectrum with the mid- 

I infrared camera ISOCAM over the wavelength range 5-16 /im and a wide range of 

Qh' orientations relative to the Sun (solar elongations 68°-113°) and the ecliptic (plane to 

O ' pole). The temperature in the ecliptic ranged from 269 K at solar elongation 68° to 244 

' K at 113°, and the polar temperature, characteristic of dust 1 AU from the Sun, is 274 K. 

■ The observed temperature is exactly as expected for large (> lO^m radius), low-albedo 

(< 0.08), rapidly-rotating, grey particles 1 AU from the Sun. Smaller particles (< 10 /xm 
^ , radius) radiate inefficiently in the infrared and are warmer than observed. We present 

H ' theoretical models for a wide range of particle size distributions and compositions; it is 

evident that the zodiacal light is produced by particles in the 10-100 fim radius range. 
In addition to the continuum, we detect a weak excess in the 9-11 fj,m range, with an 
amplitude of 6% of the continuum. The shape of the feature can be matched by a 
mixture of silicates: amorphous forsterite/olivine provides most of the continuum and 
some of the 9-11 fim silicate feature, dirty crystalline olivine provides the red wing of 
the silicate feature (and a bump at 11.35 /^m), and a hydrous silicate (montmorillonite) 
provides the blue wing of the silicate feature. The presence of hydrous silicate suggests 
the parent bodies of those particles were formed in the inner solar nebula. Large particles 
dominate the size distribution, but at least some small particles (radii ~ 1 fim) are 
required to produce the silicate emission feature. The strength of the feature may vary 
spatially, with the strongest features being at the lowest solar elongations as well as at 
high ecliptic latitudes; if confirmed, this would imply that the dust properties change 
such that dust further from the Sun has a weaker silicate feature. 
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To compare the properties of zodiacal dust to dust around other main sequence stars, 
we reanalyzed the exozodiacal light spectrum for f3 Pic to derive the shape of its silicate 
feature. The zodiacal and exozodiacal spectra are very different. The exozodiacal 
spectra are dominated by cold dust, with emission peaking in the far-infrared, while 
the zodiacal spectrum peaks around 20 fxm. We removed the debris disk continuum 
from the spectra by fitting a blackbody with a different temperature for each aperture 
(ranging from 3.7" to 27"); the resulting silicate spectra for /? Pic are identical for all 
apertures, indicating that the silicate feature arises close to the star. The shape of the 
silicate feature from /3 Pic is nearly identical to that derived from the ISO spectrum 
of 51 Oph; both exozodiacal features are very different from that of the zodiacal light. 
The exozodiacal features are roughly triangular, peaking at 10.3 /xm, while the zodiacal 
feature is more boxy. 

Keywords: zodiacal light, infrared observations, interplanetary dust 

1. Introduction 

The spectrum of the zodiacal light is a key for understanding the origin of Solar System dust 
disk and its connection to dust disks around other planetary or planet-forming stellar systems. 
Theoretically, we expect a distinction between spectra of disks around older and younger disks, as 
the originally small interstellar dust particles combine to form planetesimals, comets, and asteroids 
(Beckwith, Henning, and Nakagawa 2000). Young stars have disks dominated by small particles, 
which produce a silicate feature diagnostic of particle chemistry. Toward the well-studied star P Pic, 
the 10 /xm silicate feature was detected (Telesco et al. 1991) and its substructure measured (Knacke 
et al. 1993). Over the entire infrared spectrum of the well-observed debris disk around the star HD 
14257, several silicate features were detected, and these features are by far the dominant structures 
in the spectra (Malfait et al. 1999). For a sample of AeBe stars, silicates are generally present 
but with significant variations form star to star (Meeus et al. 2001). The asteroids and comets of 
a mature planetary system produce dust via ice sublimation and mutual collisions. While small 
particles are ejected from the solar system by radiation pressure, larger particles must be continually 
replenished to balance losses due to Poynting-Robertson drag and fragmentation due to mutual 
collisions (Griin et al. 1985). The differences in particles from star to star may be reflected in the 
variety of glass with embedded metal and sulfide (GEMS) inclusions in meteorites, some of which 
look like astronomical silicates characteristic of the diffuse interstellar medium, while others are 
more crystalline, as seen in the spectra of disks and new comets (Bradley et al. 1999). Crystalline 
silicates like SiC are present in meteorites, and they are thought to originate in the solar system since 
crystalline silicates are not observed in the ISM, possibly destroyed by supernovac shocks during 
transport from the envelopes of evolved stars (Demyk et al. 2001). Molster, Bradley, and Sitko 
(2002) suggest that large crystalline silicate grains containing iron (Fe) may have been amorphous 
grains annealed close to the sun, then transported outwards by turbulence in the early solar disk. 
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while some Fe-poor crystalline silicates in small grains may have condensed from the gas phase 
in the inner disk, then transported by radiation pressure over the solar disk to distances beyond 
Jupiter. Interplanetary dust particles captured in the Earth's atmosphere are predominantly silicate 
minerals with prominent 10 //m features in their absorption spectra (Sandford and Walker 1985). 
The spectrum of the Solar System dust is a 'laboratory' measurement, because we actually know 
some of the properties of interplanetary dust by in situ observation. The zodiacal spectrum is 
potentially a key to discriminating between old and young disks around other stars. 

The particles that produce the zodiacal light originate from both asteroids and comets. The 
mix of cometary and asteroidal dust probably varies, depending on particle size, composition, and 
orbital distribution. New comets produce crystalline particles with spectra similar to young proto- 
planetary dust disks (Crovisier et al. 1997; Hayward, Hanner, and Sekanina 2000). Analyzing 
the dynamics of particles produced by the 'mature' short-period comet 2P/Encke revealed that it 
produces only mm-sized and larger particles (Reach et al. 2000) . The broad-band infrared spectrum 
of the zodiacal light over 12-100 ^um is matched by a model with particles composed of silicates 
and the same size distribution measured in situ by Earth-orbiting satellites (Griin et al. 1985; 
Reach 1988). In addition to size and composition, orbits may distinguish asteroidal and cometary 
particles. As asteroids are relatively more confined to the ecliptic than comets, we might expect a 
difference between the mineralogy of the dust between observations in the ecliptic plane and toward 
the ecliptic pole. 

The goals of this paper are to determine whether there is a silicate feature in the zodiacal 
light and to determine whether there are spatial variations in the spectrum of the zodiacal light. 
Previous work indicated there may be some structure in the mid-infrared zodiacal light spectrum, 
but at a level of no more than 15% of the continuum. This result was obtained using the spectrum 
of the zodiacal light from 5-16 ^m measured toward one line of sight as part of the performance 
verification phase of the Infrared Space Observatory (Reach et al. 1996). The results showed a hint 
of a spectral feature with shape similar to that seen from comets (Hanner, Lynch, and Russell 1994; 
Crovisier et al. 1997; Griin et al. 2001) and the (3 Pic debris disk (Knacke et al. 1993; Fajardo-Acosta 
and Knacke 1995). Some evidence for a silicate feature in the zodiacal light was also found from 
observations with the Infrared Telescope in Space (IRTS) (Ootsubo et al. 1998: Ootsubo 2002). It 
is important not only to determine whether there is a silicate spectral feature but also to constrain 
its shape. Cometary and (3 Pic dust have already been shown to have a 10 fim. silicate feature with 
shape different from that of interstellar dust (Draine and Lee 1984) from which Solar System solids 
formed. Previous observational results for the zodiacal light were inconclusive because of potential 
detector and optical defects that had not yet been characterized, or because the spectral coverage 
did not extend beyond the red edge of the silicate feature (IRTS). In this paper, we present follow- 
up ISO observations and analysis of deep archival data that have the optimal sensitivity to the 
diffuse emission. The new observations include the ecliptic pole, ecliptic plane, and mid-Iatitudcs, 
and we compare the spectra to determine whether there is evidence for different types of particles 
overhead as compared to particles in the ecliptic, both near the Earth's orbit and exterior to it. 
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Then we compare the zodiacal light observations to exozodiacal spectra, to search for similarities 
and differences between Solar System dust and dust around other main sequence stars. 

2. Observations 

2.1. Zodiacal light: Observing strategy 

The new observations were made using the ISOCAM (Cesarsky et al. 1996) instrument on 
the Infrared Space Observatory (Kessler et al. 1996). They consist of 3 observing sequences with 
the circular- variable filter (CVF); two observations were in the ecliptic plane and one was toward 
the ecliptic pole. The observations in the ecliptic plane were performed as two-epoch, fixed-time 
observations of the same piece of sky, to measure the change of the brightness with the changing 
position of the Sun. These observations were scheduled at different times of year such that the solar 
elongations (ei) on the two dates are as close as possible to supplemental angles (ei + 62 = 180°). 
The widest separation of solar elongations makes it more possible to separate the contributions from 
the inner and outer Solar System. Taking the visibility constraints of ISO, allowing some margin 
to make it possible to schedule the observations, and choosing the ecliptic plane observations to 
be taken as far from the galactic plane as possible, the observed lines of sight were toward ecliptic 
longitude 50.1° and latitude 0.3° with solar elongations ei = 68° and €2 = 113°. The same 
astronomical observing template (AOT) configuration was used for each spectrum: 10 frames of 5 
sec duration were taken at each step of the long- wavelength circular variable filter (CVF) . We used 
the 12" pixel-field-of-view lens in order to maximize the illumination on the detector. Transient 
gain response is an important limitation of background observations like these, and the timescale 
for transient response is inversely proportional to the flux on the detector. By using the 12" pixel- 
field-of-view lens and 5 sec integration time, the flux is increased by a factor of 10 relative to 
the previous ISOCAM zodiacal light spectrum, which allowed us to reach the fainter brightness 
levels at the ecliptic pole. After each CVF observation, an observation of the same field with three 
filters (LW6, LW7, and LW8) was concatenated. These observations were used to confirm that the 
detectors were stabilized during the CVF observations. The observing conditions are summarized 
in 1. 

The ISO Short Wavelength Spectrometer (de Graauw et al. 1996) was also used in calibration 
and discretionary time to obtain deep A0T6 observations (A/ A A ^ 3000, multiple scans at the 
longest integration times) for the purpose of detecting the thermal continuum and possible spectral 
features from zodiacal dust grains. The observations were obtained in revolutions 461, 769, and 836 
at positions which corresponded to solar elongations of maximum viewable zodiacal brightnesses. 
Emission was detected in the 7 to 12 ^m range, to within the 1-a detection threshold for the SWS 
Band 2 detectors of ~50 MJy sr^^, but the sensitivity was insufhcient to allow a search for spectral 
features. No corrections to the SWS relative spectral response functions for background emission 
were ever applied, or needed as an outcome of these measurements. The observed brightness was 
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Table 1: ISOCAM Observing Log 
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consistent with that predicted by the COBE/DIRBE zodiacal light model (Kelsall et al. 1998) 
described in more detail in S2.3.1. 



2.2. Zodiacal light: Archival data 

To supplement the dedicated observations described above, we searched for archival ISO data 
suitable for deriving zodiacal light spectra. We concentrated only on ISOCAM CVF data, which 
is by far the most sensitive data for this task. We attempted to derive a zodiacal light spectrum 
from the deepest spectrometer (SWS) data, but the spectra were inadequate because the zodi- 
acal light was too faint in comparison to the dark current. One ISOCAM program in particular, 
DLUTZ.ZZULIRG, was ideal for use because (1) their intended targets were faint (< 70 mJy); (2) 
the targets were expected to be point-like; (3) the observations are very deep (over 2 hours per 
pointing); and (4) the spectra were taken aller-retour, i.e. rotating the CVF toward increasing 
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wavelengths and then back toward decreasing wavelengths. The scientific results on the intended 

targets were presented by Tran et al. (2001). Table 1 shows the observing parameters for these 
observations; we will refer to specific archival spectra by the sequential numbers listed after 'LUTZ' 
in this table. 



2.3. Zodiacal light: Data reduction 

2.3.1. Dark current and zero-point 

The observations were analyzed using routines developed at the Institut d'Astrophysique Spa- 
tiale (IAS) as part of the pre-launch calibration of ISOCAM and the performance verification (PV) 
phase of the ISO mission. For each image, we used a dark-current image obtained during the PV 
phase. To correct for the drift of the dark current between PV phase and the epoch of observation, a 
correction was obtained using the portion of ISOCAM that is not illuminated with the 12" lens and 
Fabry mirror combination. We avoided the region near the outer edge of the Fabry mirror, which 
was contaminated by stray light. The uncontaminated 'skirt' region was correlated pixel-for-pixel 
with the library dark. The library dark was offset and scaled using the zero-point and slope of this 
fit, and then it was subtracted from each image. 

The brightness at the shortest wavelengths was then then compared to the predicted brightness 
of the zodiacal light. Our model is essentially identical to the one derived by the Cosmic Back- 
ground Explorer Diffuse Infrared Background Experiment (COBE/DIRBE) (Kelsall et al. 1998). 
This three-dimensional structure of the interplanetary dust complex includes a smooth, tilted dust 
cloud, which dominates the brightness, as well as anisotropies due to the Earth's resonant dust ring 
(Dermott et al. 1994; Reach et al. 1995) and the asteroidal dust bands (Dermott et al. 1984; Reach 
et al. 1997). The wavelength-dependence is through the Planck-function kernel of the brightness 
integral as well as a set of albedoes and emissivities for each of the 10 DIRBE wavclebands spanning 
1.25-240 //m. We have generalized the model so that it predicts the brightness at an arbitrary wave- 
length (Reach 2000). We replaced these values with smooth curves for the wavelength-dependent 
emissivity and albedo. Similarly, the scattering phase-function, which was parameterized separately 
for each DIRBE waveband, was replaced by an interpolation between the phase functions DIRBE- 
derived phase functions. At ISOCAM wavelengths, scattering is only a very minor contributor to 
the brightness. Zodiacal model predictions were made for the appropriate celestial coordinates, 
observing date, and wavelength of each CVF observation. A small offset was then applied to the 
ISOCAM data, to make the brightness at 5.5-6 /xm match the COBE/DIRBE model. 



- 7- 



2.3.2. Transient response of the detectors 

These observations are at low light levels, where the observed brightness can be affected by 
transient gain response and memory effects in the detectors. We carefully checked that the detectors 
were measuring the actual spectrum of the sky as the CVF was rotated (as opposed to long time- 
constant transient gain responses). To this end, the IAS transient correction method (Coulais et 
al. 2000) was applied. This method is a straightforward calculation that does not involve fitting to 
the data and is specified by few parameters. The transient correction was verified by comparing the 
spectra taken with the CVF rotating in both the increasing- and decreasing-wavelength directions. 
The transient effect makes the detector output rise more slowly than the signal input; thus, when 
the brightness was increasing, the observed signal is lower than the true signal. Figure 1 shows the 
results of the transient correction. The first half of the observation was taken in the decreasing- 
wavelength, dccrcasing-brightncss direction, where the transient corrections were smaller. After 
correction, the decreasing- and increasing-wavelength spectra are consistent to within 1-2% of 
their average, with the decreasing-wavelength spectrum consistently brighter. This suggests that 
part of the transient gain response has not been corrected, but there is relatively little structure 
as a function of wavelength that could affect detection of spectral features. We incorporate the 
differences between the spectra taken in each direction into the systematic errors. 
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Fig. 1. — Illustration of the correction for ISOCAM CVF transient gain response, (a) The top 
panel shows the wavelength as a function of frame number for the observation of the north ecliptic 
pole. The CVF was cycled from long to short wavelength, then back up. (b) The middle panel 
shows the brightness (in electrons per second) averaged over the central part of the each frame. 
One set of dots shows the original data, and another set of dots shows the transient-corrected 
data. The detector response for decreasing signals is faster than that for increasing signals, so the 
observed and corrected data are very similar for the first half of the observation. The time constant 
for transient response is longest for low brightness levels, so the largest deviations are near the 
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beginning of the second half of the observation, where the sigla was weak and increasing, (c) The 
bottom panel shows the transient correction factor (corrected divided by observed brightness) as a 
function of frame number. 

2.3.3. Stray light removal 

A new and significant correction to the spectra was to remove stray light (Okumura 2001). This 
effect was a significant limitation of our previous attempt to measure the absolute sky spectrum 
(Reach et al. 1996). Absolute measurements are most prone to stray light because it is 
exceptionally difficult to differentiate diffuse stray light from the true sky background. Most 
flat-field calibration observations include stray light, which is desired because the stray light 
represents a real background to the detector. For the present problem of separating stray light 
from the true background, we used dedicated calibration observations of a point source observed 
over the range of CVF wavelengths and at two different locations on the array. Comparing the 
background observed through all of the ISOCAM filters, spanning 4.5-15 /xm, to the brightness 
observed with the CVF reveals systematic differences that we attribute to stray light in the CVF. 

A point source observed through the CVF generates 'ghost images' on the array; the ghosts move 
as the source is moved. The ghost images are due to multiple reflections in the optics, and these 
reflections will affect the diffiisc emission in a manner similar to the point sources. The diffuse 
stray light amplitude should be equal to the average amplitude of the ghost images for point 
sources located all over the array. In practice, we only have calibration data for two pointings of 
the point source on the array. Thus, we are forced to assume that the wavelength-dependence of 
the stray light for uniform illumination is the same as that of the average of the ghost amplitudes 
for two pointings on the array. 

For both pointings of the calibration point source, the ghost fluxes were measured and normalized 
to the main point source flux. This clearly showed the spectral dependence of the ghost intensity. 
We used the average of the ghost flux spectra for the two pointings of the calibration star to 
extrapolate the ghost intensity of the point source to the zodiacal light photometric region of the 
detector; see Blommaert et al. (2001) for further detail. ISOCAM was calibrated using the central 
point-spread-function of stars (not including the ghost images) so that extended emission (which 
includes stray light) must be corrected. The stray light amplitude depends on wavelength: it 
decreases monotonically from 35% at 6 /xm to 20% at 9.3 /xm near the edge of the 
short-wavelength CVF; then for the long-wavelength CVF it jumps to 30% at 9.5 //m, decreasing 
to 6% at 15.3 fim. The stray-light correction is discontinuous across the boundaries between the 
two CVFs at 9.5 fim. This agrees with each of our zodiacal light observations, which showed a 
significant discontinuity before stray-light correction. The stray-light correction was validated by 
comparing to data taken through 9 different ISOCAM filters. The filter data have very coarse 
spectral resolution (A/AA ~ 5), but they do not suffer as much from stray light because the filters 
are tilted with respect to the detector. To compare the filter and CVF data, we first divided each 
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intensity measurement by the DIRBE zodiacal light model (§2.3.1) appropriate for the observing 
date, direction, and wavelength. Figure 2 shows the comparison between the CVF and filter data 
(Blommaert et al. 2001). The spectra will be discussed in detail in the remainder of this paper, 
but this figure is shown in advance as a calibration validation. Before correction, there were clear 
discrepancies between the CVF and filters. After correction, the two types of observation agree 
well. Both the filters and CVF seem to show systematic deviations, at the 10% level, from the 
DIRBE zodiacal light model. 
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Fig. 2. — Comparison between the sky brightness observed through the broad, tilted filters (trian- 
gles with error bars) and through the circular- variable filter (filled circles) . Each data point presents 
the observed sky brightness after all corrections described in the text, divided by the COBE/DIRBE 
zodiacal light model appropriate for its wavelength and observing date and direction. 
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2.3.4- Spectra 

The stray-light-corrected, transient-corrected, dark-subtracted data cubes were averaged at each 
wavelength step of the CVF. The pixel-to-pixel vignetting variation was normalized by dividing 
each image by a carefully-reduced zodiacal light observation obtained during a dedicated 
calibration orbit during the mission (Biviano et al. 1998). The final spectra were generated by 
averaging over the central 12 x 12 pixels (2.4' x 2.4'). Figure 3 shows the results. The spectra are 
generally similar to each other, modulo a scale factor. They all show a similar structure in the 
9-11 /Ltm range. 



T 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 




6 8 10 12 14 16 

wavelength {/xm) 

Fig. 3. — Spectra of the zodiacal light taken on three lines of sight representing the extremes of the 

ISO operational pointing zone: in the ecliptic plane at 60° and 120° solar elongation, and toward 
the north ecliptic pole. The spectrum in the ecliptic toward e = 120° was scaled by a factor of 2.6, 
and the spectrum of the ecliptic pole was scaled by a factor of 5.5 in order to fit on the same scale 



- 12 - 



as the e = 60° spectrum. 

2. 3. 5. Reduction of archival data 

The reduction of the archival spectra used the final, revolution- and position-in-orbit-dependent 

dark current calibration for the 6" lens Blommaert et al. (2001), as well as the final spectral 
calibration and stray light correction described above. The fields were mostly clean of interstellar 
cirrus, and they spanned a useful range of ecliptic coordinates and solar elongation. Three of the 
lines of sight (numbers 16, 15, and 12 in Table 1) pass through significant amounts of interstellar 
medium. One of them (number 16) is at a relatively low galactic latitude of 13°, so interstellar 
contamination is not surprising. The other lines of sight are at high galactic latitude; however, 
two of the lines of sight (numbers 15 and 12 in Table 1) pass through high-latitude molecular 
clouds evident in far-infrared excess (Reach, Wall, and Odcgard 1998). In the analysis below, we 
will separate the three lines of sight contaminated by the interstellar medium (ISM) form the 
clean, zodiacal spectra. The spectra are shown in 4. 



-13- 




wavelength {fim) 



Fig. 4. — Spectra obtained from the deep archival observations with the 6" pixel field of view lens. 
Each spectrum is toward a different solar elongation and latitude. They were sorted in order of 
brightness and offset vertically from each other by 3 MJy sr"-*^ for clarity. 
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2.4. Exozodiacal light 

2.4.1. CAM Observations 

Wc reanalyzed the archival observations of P Pic that were first presented by Lagagc et al. (1999). 
The data were taken in a manner similar to the zodiacal light observations described above, 
except the highest magnification (1.5" pixel field of view) lens was used. The analysis of the data 
cube (images at each wavelength) was also similar. To derive the spectrum of j3 Pic, we measured 
the flux within a circular aperture (10.5" diameter) after subtracting an annular average 
background (to 22.5" diameter) at each wavelength. 

2.4.2. SWS Observations 

Observations of /? Pic have been obtained with the ISO SWS (de Graauw et al. 1996) in A0T6 

mode in programs of CWAELKEN (CRYSTAL and CRYSTAL2) and DBOCKELE (VEGAJCE), 
variously covering the full 2.4—45.2 jiva. range available between the two grating sections. Pantin 
(Waelkens) and Malfait et al. (1999) have presented the observations from the CWAELKEN 
programs. We have reprocessed the three sets of observations with the latest pipeline software 
(compatible with the ISO pipeline OLPlO.l) and interactive techniques to detect and correct the 
effects of cosmic radiation on the scan data and dark currents. These effects are manifest as 
spurious jumps in signal levels, and drastic degradation in S/N ratios over SWS bands 2 (7—12 
^m), 3 (12-39.5 /im), and 4 (39.5-45.2 /im) where the (3 Pic fiux levels arc already in the 
dark-current dominated signal regime. We have interactively treated these effects on a detector 
and scan basis, significantly improving the S/N ratios and removing artifacts. (Some such 
artifacts were visible in the earlier reductions, including band 4, where Pantin (Waelkens) and 
Malfait et al. (1999) substituted SWS Fabry-Perot data provided serendipitously but with 
unverifiable calibrations.) Then wc coaddcd the corrected products and reduced these to final 1-D 
spectra rebinned to the resolution of the instrument. This paper is focused on the ISOCAM CVF 
wavelength range. We present the 5-16 /im range of the SWS spectrum in Figure 5. We mention 
in passing the presence of CO and H2O gas absorption lines near 4.0 and 6.2 /xm, respectively, 
with strengths and amplitudes comparable to those identified in 51 Oph by van den Ancker et al. 
(2001). 

2.4.3. Spectrum j3 Pic 

Figure 5 shows the spectrum of (3 Pic from 5 to 16 /xm. The broad-band photometric points 

obtained with ISOPHOT (Heinrichsen et al. 2000) and IRAS are overlaid; they show that the 
absolute calibration of the SWS and PHOT spectra agrees. The IRTF observations of Knacke et 
al. (1993) and the ISOCAM observations are clearly different from the SWS observations, with 
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the former being significantly fainter at wavelengths beyond 10 /Ltm. The Knacke et al. (1993) 
observations cover enough of the short- wavelength, pure photosphere part of the spectrum 
(beyond the left-hand-edge of Fig. 5) that we can sec that the SWS and IRTF data agree well for 
the photosphere. The CAM data also agree in calibration with SWS at the shortest wavelength 
5.92 /xm. The differences between the spectra, which begin around 8 /xm and become very large at 
longer wavelengths, are due to the circumstellar material. The dust disk begins to dominate over 
the stellar photosphere beyond 8 /xm wavelength; by 10 /im, the circumstellar material is as bright 
as the photosphere. The aperture of the SWS and PHOT observations is significantly larger than 
that of the CAM and IRTF observations, and the colder part of the debris disk is clearly 
extended. The disk is clearly resolved in ground-based 10 /xm observations (Pantin, Lagage, and 
Artymowicz 1997), but its extent is only a few arcsec. The warm, inner disk is completely 
contained within all of the apertures discussed here. But at longer wavelengths, the disk is clearly 
extended, with a full width at half-max;imum of 8.7 x 2.9" at 25 /xm (Heinrichsen et al. 2000). 
The disk was also resolved in the longer- wavelength CAM CVF images (Lagage et al. 1999). Thus 
the CAM and IRTF observations detect primarily the inner disk, while the SWS and PHOT 
observations contain essentially the entire disk including the colder emission. 
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Fig. 5. — (a) Spectrum of f3 Pic obtained with the ISO SWS (jagged curve), CAM {X symbols) 
and IRTF (I symbols). The broadband photometry from ISO PHOT (diamonds with large error 
bars) and IRAS (asterisks with error bars) are overplotted in grey. The Kurucz model stellar 
photosphere, normalized to the 2.4-4 ^um ISO SWS flux, is shown as a solid line. Approximate 
models for the debris disk spectra are shown for the SWS spectrum (dashed line) and the CAM 
and IRTF spectra (dotted line). The differences between the observed spectra are primarily due 
to different observing apertures, as the SWS and PHOT data were taken with a large aperture 
that includes the entire cold debris disk, (h) Circumstellar spectra of the [3 Pic disk obtained by 
subtracting the stellar photosphere from the observed spectra. The curves and symbols are the 
same as in panel (a). The double-peaks at 6.4-6.8 /xm in the SWS spectrum are due to H2O. 

The lower panel in Figure 5 shows the disk spectrum, obtained by subtracting the Kurucz 1993 
ATLAS model, kp00_8250(g45) scaled to match the flux at 2.5-3 //m, for the stellar photosphere. 
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The disk spectrum is complicated, being due to a combination of particles both far from the star, 
which produces the extremely bright far-infrared emission peaking around 60 /xm (Heinrichsen et 
al. 2000), and warmer dust that produces the mid-infrared emission (Pantin, Lagage, and 
Artymowicz 1997). We will discuss these spectra and their decompositions below in § 6. 



3. Results: Continuum plus silicate feature in the zodiacal light 

The mid-infrared zodiacal spectra are dominated by a smoothly increasing brightness as a 
function of wavelength, due to thermal emission from dust at ~ 270 K, with a broad shoulder at 
9-11 fim. To remove the continuum from the spectra, we used a model for the zodiacal light that 
predicts the surface brightness as a function of wavelength for each target position and observing 
date (§2.3.1). We also fitted each spectrum to a single blackbody in order to characterize the 
continuum independent of the model. Table 2 shows the temperature and optical depths of the 
continuum spectral fits. 

Table 2: Blackbody fits and color ratios 



Spectrum 


r(K) 


loV 


^sil 


ZODY_60 


268.5 ± 0.4 


2.45 ± 0.02 


5.5 ±0.1 


ZODY_120 


244.1 ± 0.6 


1.48 ± 0.02 


5.2 ±0.7 


ZODY_NEP 


274.0 ± 1.1 


0.44 ± 0.08 


7.0 ±2.0 


LUTZ'' 


269.6 ± 5.6 


0.92 ± 0.40 


6.9 ±1.0 


LUTZ^ 


267.1 ± 0.4 


0.93 ± 0.01 


8.2 ±0.7 



average ± rms from the deep Lutz spectra (excluding those contaminated by interstellar medium) 
weighted average and statistical uncertainty from deep Lutz spectra (excluding those contaminated 
by interstellar medium) 

Fig. 6. — Comparison of the ISOCAM spectra (symbols, as in Fig. 3) with the DIRBE zodiacal 
light model for the same coordinates and date (large errors bars and solid line). 




Fig. 7. — Ratio of the observed zodiacal light spectra from Fig. 3 to a blackbody continuum. The 
continuum was fit using only wavelengths outside of 9-11 /xm. 
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Figure 6 compares the observed spectra to the their corresponding DIRBE zodiacal light model 
predictions. Figure 7 shows the observed spectra divided by their corresponding DIRBE zodiacal 
light model predictions. Because the model includes only a blackbody emissivity, the ratio 
spectrum should reveal all structures in the spectrum — both narrow lines and broad bands. It is 
clear that the model very accurate predicts the spectra: it matches to within 10%, with no 
adjusting of the model to fit the new ISOCAM data. Because of the excellent agreement between 
the COBE/DIRBE model and the ISOCAM data, and because the COBE/DIRBE data span a 
much wider parameter space of spatial and wavelength coverage, we need not delve into a detailed 
model for the zodiacal brightness distribution in this paper, but instead we refer the reader to 
Kclsall et al. (1998). The value of the new results presented here is in measuring the spectrum of 
the zodiacal light at a spectral resolution sufficient to detect spectral features. 

The ratio of the observed (ISOCAM) to model spectra shows clear systematic structures. This 
structure consists of a well-defined, broad peak at 9-11 /xm, superposed on an overall 
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'bowl'-shaped residual. The 9-11 /xm feature is the same structure as was seen in our earher 
mid-infrared zodiacal light spectrum (Reach et al. 1996), but now we have sufficient 
signal-to-noise and confidence in the instrumental calibration to confirm a detection. The 9-11 
/xm excess has also been detected in data from the Mid-Infra re Spectrometer (MIRS) on the 
Infrared Telescope in Space (IRTS) by Ootsubo et al. (1998), although the MIRS data extend in 
wavelength only to 11.59 /xm, which was not quite long enough to see the return from the feature 
to the continuum at long wavelengths. The 'bowl' shape is due to the DIRBE model not 
accurately predicting the shape of the continuum. As discussed below, the zodiacal light 
emissivity is a superposition of emission from particles of many sizes (each with their own 
temperatures) as well as particles at many locations (also having their own temperatures). The 
DIRBE model incorporates the spatial dependence of the temperature, but it assumes a single 
blackbody at each location, which is not strictly true. Because the overall shape of the continuum 
is best studied with data covering an even wider range of wavelengths than ISOCAM, we will 
emphasize the 9-11 fim feature in this paper and briefly discuss the continuum shape in 
comparing with the wider DIRBE wavelength coverage. 
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Fig. 8. — Ratio of the archival zodiacal light spectra, taken with the 6" lens and CVF of ISOCAM, 
to a zodiacal light model. The top panel shows the average of the observed to continuum ratios 
for the 14 spectra with no evidence for interstellar contamination. The bottom panel shows the 
average for the 3 spectra with interstellar contamination. The ISM lines, attributed to polycyclic 
aromatic hydrocarbons (PAH), are labeled with vertical arrows. 

Figure 8 shows the average of the archival spectra divided by their continuum from the DIRBE 
zodiacal light model. The three ISM-contaminated spectra were treated separately. They clearly 
show bright infrared bands at 6.2, 7.7, and 11.3 /xm. These are the brightest emission bands in 
spectra of diffuse interstellar matter, and they are widely attributed to polycyclic aromatic 
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hydrocarbons (PAH) (Boulanger et al. 1996). The uncontaminated spectra show no such 
interstellar lines, but instead they are dominated by the same broad band at 9-11 ^m that was 
seen in the dedicated zodiacal light observations. One straightforward result of the ISOCAM 
observations is that there are no PAH bands from the zodiacal light. The PAH bands are bright 
from the interstellar medium because of the presence of very small grains and a pervasive 
UV-bright radiation field that transiently heats the grains to high temperatures. In the Solar 
System, the interval between successive photons is shorter, and the particles will be close to their 
equilibrium temperature. This temperature is still hot enough so that the mid-infrared bands 
would be excited. The upper limit to PAH bands from the zodiacal light is < 2% of the 
continuum brightness. 

All of the zodiacal spectra, including the dedicated observations toward the ecliptic plane and 
pole (Fig. 7), the deep, archival observations (Fig. 8) and observations through the wide-band 
filters (Fig. 2), show an excess in the 9-11 /Ltm range that is characteristic of silicates. To 
characterize the brightness of the silicate feature, we define the quantity 



which is the percentage excess of the observed spectrum Ii, relative to the blackbody fit tB^{T) 
(where r and T are the fitted optical depth and temperature and B„ is the Planck function) 
averaged over the wavelengths of the feature. Table 2 shows the values of A^n for the zodiacal 
spectra. To first order, the continuum and 9-11 /im feature shapes and amplitude are comparable 
for all of the spectra. The amplitude of the 9-11 //m excess seen in the IRTS/MIRS data 
(Ootsubo 2002) is similar to that of our data. We present a more detailed comparison of the 
spectral shapes and systematic dependence on the viewing geometry in section 5.2. 



We calculated the infrared emission spectra for a range of particle compositions and size 
distributions, to learn what kinds of particles can create the observed zodiacal light spectra. 
Table 3 describes the particle compositions and optical constant data that we used. The 
calculations are an update of those that were performed to interpret the 12-100 fim broad-band 
colors of the zodiacal light (Reach 1988). We are restricting our calculations to plausible 
materials. The earlier calculations showed that pure, crystalline silicates are too cold (because 
they are nearly transparent to sunlight) while carbonaceous dust is too hot (because it absorbs 
sunlight so well). Furthermore, it is unlikely that extremely pure solids will exist in comets. 
Solids in the Solar System formed from interstellar dust, so we consider interstellar silicate dust 
as an indicator of the purity of the building blocks. A model for the optical constants of 
'astronomical silicates' (Draine and Lee 1984; Laor and Draine 1993) was found to be reasonably 




4. Comparison to theoretical models 



4.1. Compositions 
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approximated by a volume mixture of pure andesite and 3% graphite (Reach 1988) . Extensive 
infrared observations of C/Hale-Bopp demonstrated that even the sihcate dust from a new comet 
is impure: pure silicate temperatures were too low to produce the observed continuum shape and 
silicate feature ratios (Griin et al. 2001). (We confirm that this is true for many pure silicates; 
however, some amorphous silicates at 2.82 AU from the Sun have temperatures that are in the 
range [T ~ 180 K] required to explain the Hale-Bopp spectrum without need for extra 
absorption.) For this paper, we 'dirtied' the pure silicates to the same degree as the 'astronomical 
silicate' model of Draine and Lee (1984). Specifically, we use the Maxwell-Garnet mixing rule to 
add amorphous carbon with a volume filling factor f^ar (listed in Table 3) such that the Fresnel 
transmittance averaged over the wavelength range 0.2-0.7 //m (where sunlight does most of its 
heating) is at least T^jg > 0.93, the value for astronomical silicate. 

Table 3: Particle compositions and optical data 



Name 


wavelength range 


reference 


fa 
J car 


astronomical silicate (a) 


0.001-1000 Aim 


Laor and Draine (1993) 




amorphous olivine (Jena) 


0.2-500 ixm 


Dorschner et al. (1995) 





amorphous forsterite 


0.043-61.7 /im 


Scott and Dulcy (1996) 


0.1 


amorphous pyroxene 


0.2-500 //m 


Dorschner et al. (1995) 





crystalline olivine (c) 


0.01-2980 ^m. 


Mukai and Koike (1990) 


0.05 


crystalline pyroxene (d) 


6.66-487 //m 


Jager et al. (1998) 


0.033 


crystalline enstatite | 


0.043-99 /xm 


Jager et al. (1998) 


0.58 


crystalline enstatite _L 


0.043-99 /xm 


Jager et al. (1998) 


0.58 


montmorillonite (e) 


0.023-2000 /Ltm 


Koike and Hiroshi (1990) 


0.46 


glassy carbon 


0.2-100 ^m. 


Edoh (1983) 




amorphous carbon (f) 


0.105-800 ^m. 


Preibisch (1993) 




SiC 


0.001-1000 iim 


Laor and Draine (1993) 




Mg-rich crystalline olivine 


0.2-8000 /xm 


Fabian et al. (2001) 


0.37 



4.2. Absorption efficiencies 

For each particle composition, the absorption efficiency was computed on the same wavelength 
grid on which the laboratory optical measurements were provided. The computations were 
performed for particles radii on a logarithmic grid from 0.11 to 1000 /tm (particle mass 10~^^ to 
10~^ g). The calculations used Mie theory, and they were somewhat computationally intensive 
due to the inclusion of very large particles. The grain temperatures we calculated by balancing 
the solar radiation at 1 AU by the integrated thermal emission of the grains. The grains are 
assumed to be isothermal, rapidly rotating spheres; thus, absorption was assumed to occur into 
an effective cross-section tto^, and emission out of 47ra^, where a is the particle radius. The 
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largest grains can conceivably support a small temperature gradient, meaning they would emit at 
slightly higher temperatures than we have calculated. 



where dn/da is the size distribution, Bjj is the Planck function, a is the particle radius, and Qahs 
is the absorption efficiency. In practical units En is in MJy sr"-*- AU~^ or in cgs units, 
erg s^^ Hz^^ cm^^ sr^^. For our calculations, we use the following size distributions. The 
'interplanetary' and 'lunar' size distributions are based on a combined analysis of particle 
detectors on Earth-orbiting satellites, meteors, and lunar microcraters; the only difference 
between these distributions is for the smallest particles, where the 'lunar' distribution uses the 
raw lunar microcrater data, while the 'interplanetary' distribution uses an estimate based on 
/3-meteoroid production (Griin et al. 1985). The 'interplanetary' size distribution remains the best 
estimate of the size distribution of cosmic material at 1 AU from the Sun, and it has recently 
been shown to be consistent with the sizes of particles inferred from the Long Duration Exposure 
Facility (LDEF) spaceward face (McDonnell and Gardner 1998). Three other size distributions 
are included for comparison. The 'interstellar' size distribution is a power-law, dn/da oc a~^'^, as 
inferred from the interstellar extinction curve (Mathis, Rumpl, and Nordsieck 1977; Draine and 
Lee 1984); this is not intended to represent the size distribution of actual interstellar grains in the 
Solar System — the Ulysses and Galileo experiments have clearly shown that small interstellar 
grains do not penetrate into the inner Solar System (Landgraf et al. 2000) — but rather it is 
included as a simple model for comparison. The 'halley' size distribution is based on analysis of 
the particle detector data from the P/Halley encounter (Lamy, Griin, and Perrin 1987); it may 
represent the input size distribution for cometary particles that contribute to the zodiacal dust 
cloud. Similarly, the 'banner' size distribution, based on ground-based coma observations of 
short-period comets (Hanner 1984), represents freshly-produced cometary material. 

To determine which particles actually produce the zodiacal emission, we show the integrand from 
equation 1, in Figure 9. At wavelengths of strong material resonances, small particles contribute, 
but larger particles (10-100 //m radius) dominate the continuum, especially at long wavelengths. 
The size distribution is therefore a significant determinant of the strength of spectral features. 
For size distributions with relatively more small particles, the short-wavelength emission, as well 
as the spectral feature emission, can have a significant small-particle contribution, because the 
small particles are (for at least moderately-absorbing materials) hotter than the larger particles. 



4.3. 



Size distributions 



The volume emissivity is calculated by integrating over the size distribution 




(1) 
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Fig. 9. — The fractional contribution of pyroxene grains of various sizes to the zodiacal light at a 
range of wavelengths. These data are for the amorphous pyroxene model and the interplanetary size 
distribution. For each wavelength, the emissivity integrand as a function of particle size (specifically, 
{dn/da)Bi^[T{a)]7ra'^QabsiO')o- was normalized to its peak value. At most wavelengths, particles in 
the size range 10-100 ^um produce the bulk of the emission. But at the wavelengths of the silicate 
feature, a significant proportion of the emission comes from very small particles that are not opaque 
at 10 /xm. The rather complicated dependence of the emissivity on the particle size and wavelength 
makes it impossible to use single-temperature or single-size models to fit the entire spectrum. 
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4.4. Comparison of models to observations 

Model predictions for the zodiacal spectrum are shown in Figure 10. The models can be 
characterized to first order as blackbody continua plus (in many cases) silicate features at 9-11 
Hm. Using the same fitting procedure that was used on the data, we measured the color 
temperature of the contimium over the ISOCAM wavelength range; Table 4 shows the results. 
Please note that these color temperatures are only valid over the 5-16 ^m wavelength range: 
extrapolating these fits to longer wavelengths (e.g. A > 20 /im) results in very poor predictions of 
the spectrum. The carbonaceous materials are generally too hot to be the dominant, individual 
material producing the zodiacal light. The crystalline silicates were generally too cold, but after 
they were 'dirtied' by adding carbonaceous material (as described above) some of their 
temperatures arc comparable to the observed zodiacal spectra. All of the models with the Griin 
et al. (1985) interplanetary size distribution have approximately the observed color temperature, 
because large, opaque particles dominate the emission. 

Table 4: Color temperature of model continuum" 



Size Distribution 



Material 


banner 


power 


grun 


lunar 


halley-la 


hallcy-mc 


hallcy-ks 


astronomical silicate 


242.2 


250.0 


264.2 


264.9 


237.6 


247.7 


239.5 


amorphous olivine 


236.1 


248.6 


252.1 


254.9 


229.6 


244.5 


232.9 


dirty amorphous enstatite 


351.9 


420.2 


275.1 


318.8 


322.0 


399.8 


315.4 


crystalline olivine 


280.6 


289.5 


269.3 


271.9 


270.6 


296.1 


274.1 


dirty montmorillonite 


372.9 


446.2 


273.8 


325.4 


336.2 


432.3 


328.2 


glassy carbon 


404.8 


547.4 


290.5 


373.8 


361.4 


483.4 


348.3 



" temperature of best-fitting blackbody for wavelengths 6-16 jim excluding 8-13 jim 




Fig. 10. — Model spectra of interplanetary dust particles at 1 AU from the Sun, for different 
compositions and size distributions. Each panel shows the models for one particle composition 

and 5 different size distributions (long dash: halley/McDonnell, dotted: interstellar, solid: Hanner, 
dash-dot-dot-dot: halley/Lamy, dash-dot: lunar, dashed: interplanetary; listed in order of silicate 
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feature strength for astronomical silicate). The 'power- law' size distribution, which has the largest 
abundance of small particles, has the highest-contrast spectral features, while the Griin et al. (1985) 
interplanetary size distribution, which has the lowest abundance of small particles, has essentially 
no silicate feature. Panel (a) represents interstellar material; panels (h), (c) represent amorphous 
silicates; and panel (d) \s a. crystalline silicate; panel (e) is & meteorite sample; and panel (f) is an 
amorphous carbon sample. In each panel, the average observed zodiacal spectrum is overplotted 
as a set of thick black circles. It is clear that there arc no high-contrast features in the observed 
spectrum, ruling out essentially all of the small-particle-cominated size distributions. 

The observed zodiacal spectrum has very little structure; specifically, there are no features above 
10% of the continuum. In order to assess which model spectra show comparable features, we 
divided each model spectrum by its blackbody continuum fit (see Fig. 11). First, it is evident that 
all of the spectra, both model and observed, exhibit a rise at short wavelength relative to their 
blackbody fits. This rise can be easily explained by a range of dust temperatures contributing to 
the emission. The blackbody fit will yield the color temperature that dominates the 6-15 jim. 
wavelength range. It will not include hotter dust, which dominates at shorter wavelengths 
because of the steepness of the Wien portion of the Planck curve. The short-wavelength rise can 
be due to small silicates or to carbonaceous material. 
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Fig. 11. — Model spectra from Fig. 10, divided by a blackbody fit to their mid- infrared continua. 
The fits were performed in the same way for the observed data, which are overplotted on each 
panel. The rise in the observed ratio spectrum at short wavelengths is reproduced in the model 
spectra; it is due to the spread in temperatures as a function of particle size. The bump at 9-11 
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/xm that could be a silicate feature is roughly comparable to those predicted by the silicate models, 
but it does not match any of these models in detail. The zodiacal dust must include a range of 
different mineralogies; they cannot be dominated by a single crystalline or amorphous silicate. 

The spectral features predicted by the models are generally much too strong and too structured 
to match the observed spectrum, except for the large-particlc-dominatcd models. In Figure 11, it 
is clear that the individual silicates, both crystalline and amorphous, generate distinctive shapes 
in their 9-11 /xm spectral features. The observed spectral feature is rather broad and smooth, 
lacking sharp peaks. However, the observed line is 'boxy' and the short- and long-wavelength 
edges of the observed feature roughly correspond to resonances in silicates. We can reach several 
simple conclusions from these observations. First, the short-wavelength excess emission seen in the 
observed spectrum is more likely due to carbonaceous (or another strongly-absorbing) material; if 
it were due to small, hot silicates, then the 9-11 /um feature would be much stronger. Second, the 
zodiacal light is produced mostly by large (> 10 /xm) particles, and the size distributions with 
large abundances of small particles can be ruled out. Prom Table 5, only the interplanetary and 
lunar size distributions from Griin et al. (1985) are plausible. The other size distributions, which 
are based on observations of comets and the interstellar medium, all have far too many small 
particles, making their spectral features and too prominent (i.e. with Deltagu much larger than 
observed). Thus is is clear that zodiacal dust is very different, in size or composition, from the 
dust that produces the observed comelary or interstellar silicate spectral features. 

Table 5: Particle size distributions" 



Name 


^small 


^large 


reference 


interplanetary 


-2.16 


-2.00 


Griin et al. (1985) 


lunar 


-3.74 


-2.05 


Griin et al. (1985) 


interstellar 


-3.50 


-3.50 


Mathis, Rumpl, and Nordsieck (1977) 


Halley/Lamy 


-2.53 


-3.40 


Lamy, Griin, and Perrin (1987) 


Halley/McDonneU 


-3.11 


-3.96 


McDonnell et al. (1991) 


Hanner 


-2.66 


-3.59 


Hanner (1984) 



" ksmall and harge are the approximate power-law indices of the size distributions between 
particle sizes 0.3-2 fim and 2-20 /xm, respectively. 



4.5. Mixture model 

We constructed a non-unique, but illustrative, model that quantitatively matches the broad 
zodiacal spectrum. The model uses a large-particle-dominated size distribution with just enough 
small silicates to make the weak 9-11 /xm feature. A better fit to the detailed shape of the silicate 
feature using different particle shapes is described below; the model described here addresses the 
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overall spectral shape including feature and continuum. First, because the interplanetary size 
distribution yields no silicate features and the lunar (and other) size distributions make too 
strong a silicate feature, we had to construct a new size distribution. This was simply done by 
taking the difference between the model spectra for the interplanetary and lunar distributions, 
and adding in a fraction of that difference to the interplanetary distribution until it makes silicate 
feature amplitude comparable to the observed one. This is valid because the model is a linear 
function of size distribution (eq. 1), and the lunar size distribution is identical to the 
interplanetary size distribution except for small particles (Griin et al. 1985). Second, we had to 
combine different materials in order to generate a wide silicate feature. 




6 8 10 12 14 16 

wavelength (/xm) 



-32- 



Fig. 12. — Comparison of a mix of silicates and carbonaceous materials to the observed spectrum 
of the zodiacal light. The top panel shows the observed spectrum (circles with error bars), together 
with the mixture model (solid line) and the brightness of each of the 3 components of the mixture 
(dash: dirty amorphous olivine, dash-dot: glassy carbon, dash: crystalline olivine). The bottom 
panel shows the ratio of the intensity to a blackbody continuum fit for the observed (symbols with 
error bars) and mixture model (line) spectra. 

Figure 12 shows the mixture, compared to the observed spectrum. The fit to the spectral feature 
is relatively poor, but the main features of the observed spectrum are reproduced. The mixture 
contains three materials: amorphous forsterite/olivine from Scott and Duley (1996) [80% of the 
13-15 /im intensity], crystalline olivine from Mukai and Koike (1990) [10%], and amorphous 
carbon from Edoh (1983) [10%]. Relative to cosmic elemental abundances, this is a silicate and 
oxygen-rich mixture. Most of the carbon is mixed in with the amorphous forsterite/olivine (Scott 
and Duley 1996), which was 'dirtied' using a 10% volume mixture of glassy carbon; the free, 
amorphous carbon grains contain the remaining 35% of the carbon. The relative paucity of 
carbonaceous material could be due to its volatility: the grains, whether originally cometary or 
asteroidal, have been exposed to sunlight for > 10^ yr (Griin et al. 1985), which may be long 
enough to remove their volatiles. The elemental abundances from our mixtTirc model should not 
be taken as the total elemental abundances of the interplanetary dust, because they are weighted 
toward particles that are relatively bright in the mid-infrared. The elemental abundances from 
our mixture model are different from those of fresh cometary ejecta, which has a significant 
gaseous component rich in the lighter elements (e.g. as detected by impact detectors that flew 
through comet Halley's coma in 1986). The size distribution is dominated by large (10-100 /Ltm 
radius) particles. For crystalline olivine and amorphous carbon, the size distribution was 
intermediate (50% each) between the interplanetary and lunar size distributions, while the 
amorphous olivine, which dominates the spectrum, had a size distribution very close to 
interplanetary (92%, with 8% lunar). 

In order to fit the detailed shape of the silicate feature, we had to consider the effect of particle 
shape. Inspecting Figure 11 carefully, it is evident that none of the theoretical models for spheres 
produce silicate features extending beyond 11 /xm, as the observed feature does. This is true 
despite the presence of strong resonances in crystalline olivine at 11.4 iim (Fabian et al. 2001). It 
is a well-known property of spheres that their absorption cross-sections do not necessarily peak at 
the wavelengths of the resonances, instead having wavelength dependence, in the small particle 
approximation, 

Qsphere _ . j ^ ~ ^ 
^abs ~ "i-i^I" ^ _^ 2 ' 

where x = 27ra/A, e = m^, and m is the (wavelength-dependent) complex index of refraction. This 
equation combines the real and imaginary parts of the index of refraction in a very specific, 
nonlinear way. The peaks of silicate features from spherical particles are significantly shifted from 
other shapes (Yanamandra- Fisher and Hanner 1999). To generalize the theoretical calculations, 
we fitted the shape of the silicate feature using the 'continuous distribution of ellipsoids' (CDE) 
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approximation, in which ah possible axis ratios have equal probability; the absorption efficiency 
averaged over shapes is then 



which does not shift emission features from resonances as much, though it does tend to broaden 
them (Bohren and Huffman 1983). This equation only applies in the small-particle limit, and the 
full equations for a distribution of larger ellipsoids are not known and probably computationally 
very challenging. Since the particles producing the continuum of zodiacal light are large (as 
discussed above), wc cannot at present use a full theoretical CDE calculation (solving for the 
temperature for each particle size) for the zodiacal light. More sophisticated models of 
non-spherical particles, such as the discrete dipole approximation (Draine 1988) or T-matrix 
method (Mischenko 1991) would be able to treat the problem self-consistently, but they are far 
beyond the scope of this project. However, the zodiacal light silicate feature is produced by 
particles small enough that the small-particle limit can be used. Therefore, we continue to use the 
models for spheres, as presented in the previous sections, to calculate the continuum spectrum. 
To model the silicate feature, we calculated CDE absorption efficiencies for all of the materials 
and for 0.4 /xm particle radii (see Fig. 9). A linear regression of the absorption efficiencies was 
fitted to the observed spectrum. If the various silicate minerals have the same size distribution, 
then the scale factor from the regression analysis gives the relative abundance of each mineral. 

The best fit to the silicate feature includes three components; Figure 13 shows the mixture and 
each component. This solution was found by trying each three-way combination of 40 different 
theoretical models (each of the materials described above, some with and without different 
amounts of carbonaceous material), and identifying the solution with minimum x^- (Solutions 
with negative coefficients were discarded as nonphysical, leaving 3741 models involving 
combinations of various materials.) The most abundant silicate is amorphous olivine (with a scale 
factor of 1.12), which is the best single- mineral fit to the spectrum. However, amorphous olivine 
cannot fit the blue or red wing of the observed feature, which is distinctly 'boxy.' The blue wing 
requires a different silicate mineral, and a reasonable fit can be obtained with a hydrous silicate 
(montmorillonite, with a scale factor of 0.25). There is some indication of a more detailed 
substructure in the blue wing that can be fitted with this mineral, with a possible hump in the 
spectrum around 8.6 /xm. None of the minerals we considered have this structure, but it is 
relatively subtle and uncertain so we did not pursue it further. The red wing of the zodiacal 
feature, including the small peak at 11.35 /xm, is reasonably fit by crystalline olivine (with a scale 
factor of 0.12). 
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Fig. 13. — Comparison of a mixture of silicates with all possible ellipsoidal shapes to the observed, 
continuum-subtracted, spectrum of the zodiacal light. The solid line is the mixture fit, the dotted 
line is the amorphous olivine contribution, the dashed line is the montmorillonite contribution, and 
the dash-dotted line is the crystalline olivine contribution. 

The presence of a hydrous silicate in the mixture required to fit the shape of the silicate feature 
indicates that the parent bodies of the zodiacal particles were formed in a location where water 
could be mixed with crystallizing minerals. This may suggest formation in the inner Solar System, 
where temperatures and pressures were high, or processing on an active surface (or in collisions). 
A mixture of water and minerals is certainly expected for comets, but the water is less likely to be 
contained so deeply within the silicate minerals. These subtle clues can reveal relationships 
between zodiacal particles and their parent bodies, and they need further exploration. 

Using the best-fitting model, we can predict the spectrum at wavelengths outside the region we 
observed. Due to the prevalence of large particles, and the weakness of the 9-11 /xm feature, we 
predict a relatively featureless spectrum. The CDE fit includes crystalline olivine, which has 
features at 20, 21, 24, and 35 /um, and montmorillonite, which has features at 19, 23, and 29.5 /xm 
(brightest features in italics). We note that earlier measurements of montmorillonite (Koike, 
Hasegawa, and Hattori 1982) also show a broad feature in the 45 to 55 jum range that is not 
reproduced in their later derivations of the optical constants (which we use here) from Lorentz 
oscillators. The amplitude of these features is not accurately predicted by the CDE model, 
because it is a single-size model. A rough estimate can be made by using the CDE mixture model 
to give the detailed features and the spherical model to give the continuum level. Interestingly, 
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the strongest crystalline olivine and hydrous silicate features alternate in wavelength, such that 
the sum is more like a lumpy continuum. Only the 9-11 fim. feature is strong in the mixture. The 
next brightest features (in italics above) arc only 1 % of the continuum. Thus, apart from the 
9-11 iim feature, we do not expect significant structure in the zodiacal light spectra obtained 
with the Infrared Spectrograph on the Space Infrared Telescope Facility (Houck 1995). By 
extension, we predict a relatively featureless spectrum for the terrestrial portion of debris disks 
around mature stars, where the primordial dust has been removed and the particles are derived 
from asteroid collisions and comet fragmentation. 



5. Relation between zodiacal light and parent bodies 

5.1. Comparison to IDPs and cometary dust 

Sandford (1987) summarized the properties of interplanetary dust particles (IDPs) collected in 
the atmosphere. He divided the particles into three classes: olivines, pyroxenes, and layer-lattice 

silicates. All of the IDPs arc heterogeneous, with the carbonaceous material that survived 
atmospheric entry spread throughout the grains. The collected particles are rather porous, and 
the 'empty' regions were probably at least partially filled with more volatile materials before 
atmospheric entry. These properties are consistent with our 'dirtied' silicate sphere models, 
although the actual porosity of the grains in space is not known and could affect the spectral 
properties. The infrared transmission of a sample of IDPs was measured by Sandford and Walker 
(1985). All of the silicates had 9-11 fim features. The olivine particles presented substructures at 
9.3, 10.1, and 11.2 ^m, while the pyroxenes presented substructures within the 9-11 /im band 
that varied strongly from particle to particle. The layer-lattice silicates have strong 0-H features 
that would have been observable in the zodiacal spectrum. Thus, small particles composed of 
layer-lattice silicates are unlikely to be very abundant in interplanetary space. The observed 
zodiacal light spectrum, on the other hand, bears a strong resemblance to the olivine and 
pyroxene IDPs. The strong particle-to-particle mineralogy variations would, when blended 
together along the line of sight in a zodiacal light observation, lead to a relatively broad and 
smooth feature, which is in fact what we observe. 

Comets in the inner Solar System present relatively strong 9-11 /Lim silicate features (Hanner, 
Lynch, and Russell 1994). These features vary from comet to comet, and they may be stronger for 
new comets. The best-observed cometary silicate feature is for comet Hale-Bopp, both because it 
was very bright and because of advances in infrared detection. Wooden et al. (1999) presented 
extensive ground-based observations of the 9-11 fim Hale-Bopp, and Crovisier et al. (1997) 
presented the ISO spectrum over a wider wavelength range. Comparing the continuum-subtracted 
zodiacal light and Hale-Bopp silicate features (shown together in Fig. 15 below), the amplitude of 
the silicate feature from the Hale-Bopp dust is twice the brightness of the underlying continuum. 
In our nomenclature, Agn ranged from 160 to 200 as the Hale-Bopp moved from 2.8 to 0.9 AU 
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from the Sun (Wooden et al. 1999). Comets Halley and Mueller showed similarly strong silicate 
features, while comets BorrcUy and Faye had weaker features, with A^j; ~ 20 (Hanner, Lynch, 
and Russell 1994). The comctary silicate features are several to 20 times stronger than the silicate 
feature in the zodiacal light. For comets with infrared spectra, the particles are inferred to be 
much smaller than those that dominate the zodiacal light, due both to their bright silicate 
features and their continuum temperatures being much higher than blackbodies at the 
corresponding distance from the Sun. These coma grains are so small that they feel enough 
radiation pressure to put them in hyperbolic orbits. It is clear that these small dust grains are not 
the origin of interplanetary dust, unless the grains can grow in the interplanetary medium. 
Molster, Bradley, and Sitko (2002) suggested that some Fe-poor crystalline silicates may have 
condensed out of the gas phase close to the Sun, and the small grains transported by radiation 
pressure over the solar disk to distances beyond Jupiter. Larger cometary particles have been 
detected by in situ observations at P/Halley (McDonnell et al. 1991) and inferred by dynamic 
modeling of the morphology of the dust around C/Austin (Lisse et al. 1994) and P/Encke (Reach 
et al. 2000). Let us assume that the small comctary grains have similar mineralogy to the larger 
ones, so it is useful to compare the mineralogy of the cometary and zodiacal particles. First, 
comet grains appear to be more crystalline, with sharper substructures in their 9-11 feature 
(Wooden et al. 1999), and strong crystalline features in the 35 /im region (Crovisier et al. 1997). 
In Figure 15, the 11.3 jim crystalline olivine peak is clear in the Hale-Bopp silicate feature. 
Second, the comet grains appear to be Mg-rich (Wooden et al. 1999), which is similar to our 
results for the zodiacal light, where a forsterite material (Mg-rich olivine) matches the spectrum 
better than fayalite (Fe-rich) . We cannot say with confidence whether there is a real distinction 
between what seem to be olivines (Si04) matching the zodiacal spectrum and the mixture of 
olivines and pyroxenes (SiOs) matching the Hale-Bopp spectra (Wooden et al. 1999), because the 
zodiacal spectrum is too dominated by large particles to present clearly-defined spectral features. 

5.2. Spatial variations of the zodiacal spectrum: cometary versus asteroidal origin 

for dust 

A goal of this project was to determine how the zodiacal spectrum varies from place to place. To 
first order, the spectra are remarkably similar. The brightness is clearly higher at low ecliptic 
latitudes and toward the Sun, as shown in the definitive maps of the zodiacal cloud by 
C05E/DIRBE (Kelsall et al. 1998). For a more detailed comparison. Figure 14 shows the 
temperature and silicate excess determined from fits to the ISOCAM CVF spectra described 
earlier in this paper. The fit values are displayed in the relevant Solar System coordinate system; 
specifically, they are plotted as a function of angular separation from the Sun and absolute 
ecliptic latitude. The spectra with the lowest color temperature are those taken farthest from the 
Sun and at relatively lower ecliptic latitude. The warmest spectra are those closest to the Sun 
and at high ecliptic latitude. These variations are intuitively expected, because the sight-lines 
that contain the coldest dust will be those that see farthest toward the outer Solar System. The 
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spectra taken closest to the Sun sample the hottest dust (reaching as close as 0.9 AU). The reason 

the polar spectrum is among the warmest is that it contains almost exclusively dust at 1 AU, 
while the ecliptic plane spectrum closest to the Sun samples dust at a range of distances, from 0.9 
AU outward to distances where the dust is cooler. 
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Fig. 14. — The color temperature (a) and silicate excess (h) derived from ISOCAM spectra of 
the zodiacal light. Both the archival and dedicated zodiacal light observations are included; the 
dedicated observations are distinguished by larger font. 
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The value of the polar temperature is very similar to the temperature (278 K) expected for 
rapidly-rotating, grey grains at 1 AU from the Sun. This clearly demonstrates that the spectrum 
is dominated by grains large enough that they emit efficiently in the mid-infrarcd. Specifically, 
2'Ka/X > 1, where A ~ 30 nm is the longest wavelength where the grains must emit efficiently in 
order to behave as grcybodies. This sets a limit on the grain radii of a > 5 /xm, consistent with 
the theoretical models described above, which have the emission produced by particles in the 
10-100 /um size range. 

The variation of the color temperature with viewing direction is qualitatively as predicted by the 
DIRBE zodiacal light model (Kelsall et al. 1998). To test this quantitatively, we evaluated a 
generalized version of the DIRBE zodiacal light model (§2.3.1) over the mid-infrared wavelengths 
at which the ISOCAM observations were made, and we used the same color-temperature-fitting 
routine as we used on the data. The model yields colors with approximately the same variation of 
temperatures with solar elongation and latitude. The pole-to-plane temperature difference at 
solar elongation e = 113° is 30 K from ISOCAM and 24 K from the DIRBE model; the in-plane 
(e = 66° - 113°) difference is 25 K from ISOCAM and 23 K from the DIRBE model. But the 
DIRBE model is consistently warmer than the ISOCAM spectra, by ~ 6 K. The discrepancy can 
be traced to the high dust temperature at 1 AU, Tq = 286 K, used in the DIRBE model. This 
temperature used in the DIRBE model but was not included as a free parameter; an independent, 
DIRBE-based model by Wright (1998) also adopted a value of Tq, but the adopted value was 
To = 280 K, which is in better agreement with the ISOCAM data. 

While the continuum temperature was expected to show a reflex variation as a function of 
viewing direction, due to changing dust temperature as a function of distance from the Sun, the 
amplitude of the silicate feature, expressed as a fraction of the continuum, is expected to remain 
relatively constant within the inner Solar System unless the properties of the material change 
from place to place. In Figure 14, the silicate feature amplitude shows a possible trend as a 
function of viewing direction. Specifically, the feature is somewhat stronger at lower elongations 
and higher latitudes. This suggests that dust relatively farther from the Sun has a weaker silicate 
feature than dust within ~ 1.1 AU of the Sun, and possibly also dust with higher- inclination 
orbits has a somewhat stronger silicate feature. 

The variation of the silicate feature is small, if present, which suggests that the mix of dust grains 
does not vary dramatically as a function of distance from the Sun or from the ecliptic. This 

relative uniformity applies over the range of distances within which the zodiacal light is produced. 
The radial range is not great: at 12 fim, in the ecliptic, toward 90° elongation, 90% of the 
zodiacal light is produced between the orbits of Earth and Mars (Reach 1991). 

A variation in the apparent mineralogy of dust in the coma of C/Halc-Bopp was detected by 
Wooden et al. (1999), who observed the silicate feature from the comet at a range of heliocentric 
distances. They found that both the amplitude of the feature (relative to the continuum) and its 
shape varied with heliocentric distance. They inferred that there were multiple silicate 
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components (similar to our three-component fit to the zodiacal spectrum), and that the 
temperature of at least one of the components was distinct from the others. The coldest 
component does not become bright enough to contribute to the 9-11 /um feature until it is 
relatively closer to the Sun. The net effect is that the line-to-continuum ratio of the silicate 
feature increases as the heliocentric distance decreases. This effect goes in the same sense as the 
effect seen in the zodiacal light, suggesting a similar mechanism may be at play. 

There have been other suggestions of variations of dust properties as a function of location in the 
Solar System. Renard, Levasseur-Regourd, and Dumont (1995) found that the dust albedo 
decreases as a function of distance from the Sun as ^ oc r~°-^^. That is, they found that grains 
are more dark at greater distances from the Sun. This result could be qualitatively in accord with 
a decrease in the amplitude of the silicate feature, if the grains are larger at greater distances 
from the Sun, because larger grains have less silicate feature and are relatively darker. (Larger 
grains are more efficient at infrared absorption per unit visible absorption due to physical optics.) 
Whether the albedo actually varies is extremely difficult to say with confidence, and the same 
applies to our suggested variation in the silicate feature. Prom remote sensing, we only retrieve an 
integral along the line of sight, which cannot be uniquely inverted. The results from the Helios 
spaceprobe, which measured the brightness, color, and polarization of the zodiacal light while 
orbiting between 0.3 and 1 AU from the Sun, are considerably more reliable. Leinert et al. (1981) 
showed that the Helios results require a radial variation of the volumetric scattering cross-section 
nA oc r-^-^=^° °5. The Kelsall et al. (1998) fit to the DIRBE infrared data yielded a density profile 
n oc 7.-i-34±.02 consistent with the volumetric scattering cross-section from Helios; 

combining these results, no variation of the albedo is required. However, the infrared model 
assumed constant dust properties, and it has not been determined whether a model with variable 
dust properties fits the data significantly better. The only convincing evidence for a variation in 
dust properties is the Helios observation that dust closer to the Sun generates relatively less 
polarization than dust at 1 AU (Leinert et al. 1981). 

Radial variations in dust properties, if present, can be explained theoretically. The observations 
have suggested that dust is somewhat darker, more polarizing, and less silicate-feature-producing 
at larger distances from the Sun. Greenberg and Hage (1990) predicted varying properties of 
cometary dust, with relatively larger number of small grains in the inner Solar System being due 
to disintegration of cometary particles that consist of many small grains stuck together by 
volatiles, which evaporate as the grains approach the Sun. The same qualitative result is 
produced by grain-grain collisions, regardless of particle morphology, as very large grains, 
spiraling inward from the asteroid belt (for example) under the influence of Poynting-Robertson 
drag, are shattered by collisions before they reach the inner Solar System. Griin et al. (1985) 
calculated the coUisional lifetimes to be shorter than the Poynting-Robertson lifetimes for grains 
larger than 100 //m. Asteroidal grains of this size would not survive into the inner Solar System. 

The carriers of the silicate feature are actually very small particles. Figure 9 showed that the 
feature is produced most efficiently by submicron particles. Such are predicted to have very short 
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lifetimes, because solar radiation pressure exceeds gravity (Burns, Lamy, and Soter 1979). For the 
materials used in our mixture model, we calculated the ratio of radiation pressure to gravitational 
force, /3. If the particles are compact, with volume densities p ~ 2 g cm^'^ such as typical of 
terrestrial rocks and stratospheric-collected IDPs (Love, Joswiak, and Brownlee 1994), then 
radiation pressure exceeds gravity (/? > 1) for silicate particles (both crystalline and amorphous) 
with radius smaller than ~ 0.4 /xm. (Amorphous carbon particles always have /? < 1, unless their 
densities are low.) If the particles are porous, then radiation pressure exceeds gravity for a 
proportionally larger radius, because /? oc l/{pa). Such particles should be relatively rapidly 
ejected from the Solar System. However, it is clear that particles of this size are present in the 
inner Solar System, as they were observed by dust detectors aboard Pioneers 8 and 9 and 
HEOS-2 (Griin et al. 1985). There is a continuous source of these so-called /?-meteoroids, from 
the coUisional comminution of larger particles that are spiraling toward the Sun under the 
influence of Poynting-Robertson drag. Thus both comets and asteroids particles can be sources of 
both large and small particles in the inner Solar System: Poynting-Robertson drag brings them 
in, and collisions grind them down. 

To distinguish a cometary and asteroidal origin for the zodiacal dust, a telling variation could be 
a dependence of the strength (or shape) of the silicate feature on ecliptic latitude. While asteroids 
are confined to the ecliptic, comets have a wider range. New comets arrive at essentially all 
inclinations; however, nearly all of the dust they produce leaves the Solar System due to the high 
eccentricity of the comet's orbit and the effects of radiation pressure (Burns, Lamy, and Soter 
1979). Short-period comets arc the most plausible cometary source of zodiacal dust. The mean 
inclination of main-belt asteroid orbits is 7.9°, while the mean inclination of comet orbits is 19° 
(Cox 2000). Thus dust from short-period comets would be on relatively higher-inclination orbits 
than asteroidal dust, and it would be relatively more prevalent toward higher ecliptic latitudes. 

Can we see systematic differences in the silicate feature versus ecliptic latitude? In Figure 14, 
there is no strong trend. The high-latitude spectra do show a silicate feature, and that feature 
appears to have approximately average amplitude, though perhaps on the high side. We suspect 
that any dependence on the orbital inclination of the parent body is masked by radial variations, 
so that from our vantage point near the ecliptic plane we cannot easily separate radial from 
vertical variations. The best chance for such a separation is to consider the spectra along the 
meridian where the solar elongation is 90°, and looking only at latitudes above 10°. Although the 
ecliptic plane will clearly contain relatively more distant dust, by moving along this meridian we 
will keep the mean solar distance to the emitting region at least roughly constant. We find that 
the high latitude points have an silicate excess Agu = 7 it 1, while the low-latitude points have 
Agii = 7 lb 1. Thus there is no evidence for a segregation of dust properties as a function of 
ecliptic latitude. Our result (or non-result) is still consistent with estimates of parent sources for 
the interplanetary dust, which show that both comets and asteroids arc viable sources. For 
example, Whipple (1955) predicted that short-period comets produce enough dust to maintain 
the zodiacal cloud, while Grogan, Dermott, and Durda (2001) showed that the asteroid families 
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contribute a third of the cloud and non-family asteroids can produce the rest. 

6. Comparison of zodiacal and exozodiacal spectra 

The new mid-infrared spectra obtained from ISO allow us to compare the properties of material 
in the interplanetary medium, comets, and in debris disks around other stars. Figure 15 compares 
the spectrum of the zodiacal light, comet Hale-Bopp (Crovisier et al. 1997), the (3 Pic debris disk, 
and the interstellar medium (Draine and Lee 1984). To make this comparison, we have subtracted 
smooth continua from each spectrum. The zodiacal light continuum was obtained by first 
removing the Kelsall et al. (1998) zodiacal light prediction for this line of sight and observing 
date. (This model uses a blackbody kernel so it will not affect the silicate feature shape). Then 
wc subtracted a second-order polynomial in the 5-16 /im range excluding 8-12 ^um. For 
Halc-Bopp, the continuum is simply a second-order polynomial fitted over this same wavelength 
range. For the interstellar spectrum, we used the theoretical model described above for 
astronomical silicate particles heated by the Sun at 1 AU with the interstellar size distribution 
(Mathis, Rumpl, and Nordsieck 1977), and then subtracted a second-order polynomial fitted over 
the same wavelength range as for the other spectra. 
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Fig. 15. — Continuum-subtracted spectra of (a) the zodiacal light, (b) comet Hale-Bopp (Crovisier 
et al. 1997), (c) the (3 Pic debris disk, and (d) interstellar particles (Draine and Lee 1984). The 
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spectra are roughly in evolutionary order, with the initial building blocks from the interstellar 
medium at the bottom, then the agglomerating particles in the /3 Pic disk, then the comet, then 
the cometary and asteroidal fragments of the zodiacal light. In panel (c), the observed spectra 
(same curve and symbol coding as Fig. 5) minus the photosphere and continuum debris disk are 
plotted. Comparing the zodiacal and exozodiacal features, the former appears redder and rela- 
tively boxy (like the cometary spectrum) , while the /3 Pic feature appears smoother and somewhat 
triangular. Vertical dashed lines arc shown at two suggestive wavelengths: the 11.3 /xm peak in 
the cometary spectrum (attributed to olivine) and the 9.4 /im peak of the interstellar astronomical 
silicate spectrum. (In panel (c), the emission features in the 6.4-6.8 /xm range are due to II2O.) 

For (3 Pic, the continuum was modeled using a stellar atmosphere model for the star and a simple 
approximation of the disk. The separation of the silicate feature and continuum is difficult for (3 
Pic because the feature is so smooth, and there is such a wide range of physical conditions 
contributing along the line of sight — from the outer disk, which dominates at long wavelengths 
(Heinrichsen et al. 2000), to the inner disk, which is only revealed at mid-infrared wavelengths 
(Pantin, Lagage, and Artymowicz 1997) and high-resolution, coronagraphic optical images (Heap 
et al. 2000). For the large-aperture observations, which include the entire debris disk, we made a 
simultaneous fit to the continuum in the SWS and PHOT data (5-200 /xm) using two modified 
blackbody components: 

= [(l.OTzbO.ll) X 10-^5^(77.5 ±1.4 K) (2.37 ±0.18 x 10-^)5^(197.5 ± 2.5 K)](12/Ltm/A) (2) 

where Bi,{T) is the Planck function. The two components represent the 'outer' and 'inner' 
exozodiacal light, respectively, with the 'outer' component being the one that dominates the 
far-infrared brightness and the mass of the disk. This fit, shown as a dashed line in Figure 5, 
becomes brighter than the photosphere at wavelengths longer than 12 ^um. For the smaller 
aperture observations by ISOCAM, we approximated the continuum with a single blackbody 
component (hence the fit is only valid at wavelengths shorter than 15 /xm), 

= 1.54 X 10-^5^,(300 K), (3) 

which is shown as a dotted line in Figure 5. The continuum is weaker in the ISOCAM and 
ground-based mid-infrared spectra, compared to the SWS and PHOT data, because part of the 
'outer' disk is resolved and excluded from the former. The comparable temperatures of the 
ISOCAM continuum and the 'inner' component of the SWS/PHOT continuum fit suggest they 
both arise from similar parts of the disk, but even the 'inner' SWS/PHOT continuum is produced 
from a somewhat greater distance from the star than the CAM/ground-based continuum. Recent 
high-angular-resolution observations with Keck by Weinberger et al. (2002) resolve the thermal 
emission from the disk even better than those described here, revealing that the spectral shape 
changes with distance from the star, such that the silicate feature is largely produced from the 
inner 1" of the disk. A more detailed model for the debris disk is needed to study the shape of the 
continuum feature in detail; for our purposes, the simple models are sufficient to allow us to 
extract the silicate feature. 
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The (3 Pic silicate feature appears similar to that found in recently-published exozodiacal light 
studies. After the first detection of the silicate feature from (3 Pic, Telesco et al. (1991) had 
already noted that the silicate feature is wider than that the interstellar feature (Draine and Lee 
1984). Waelkens et al. (1996) showed that the silicate features of (3 Pic and 51 Oph are very 
similar, van den Ancker et al. (2001) presented the remarkable spectrum of 51 Oph, for which the 
debris disk is so bright that subtraction of the central star is hardly needed in the mid-infrared. 
They noted the extremely wide, and somewhat triangular, shape of the silicate feature, and 
hypothesized that the red wing of the feature, which is wider than silicate minerals we have 
studied to date, may be due to a blend with a second mineral. They also found that the 51 Oph 
silicate feature is identical to that of symbiotic star V385 Cen but wider (on the red side) than 
that of Herbig Ae star AB Aur, leading to the tantalizing speculation that the dust in 51 Oph 
may be " 'fresh', i.e. recently formed rather than accreted from the interstellar medium and 
processed in a proto-planetary disk" (van den Ancker et al. 2001). Given the extreme similarity of 
silicate features between [3 Pic and 51 Oph, it is likely that the broad, triangular shape may be a 
characteristic feature of proto-planctary disks, and may in fact be due to processing (including 
agglomeration) occurring in the disk. 

It is clear that the overall spectrum of the zodiacal light is very different from exozodiacal light as 
evidenced by (3 Pic or 51 Oph. These exozodiacal spectra are dominated by dust much further 
from their exciting stars and therefore much colder. The 10 /xm silicate feature from of these disks 
is produced by dust in the inner disks. Thus the line-to-continuum ratio (which is a good 
indicator of particle size and mineralogy) is very difficult to compare: one would have to isolate 
the contributions form the outer and inner disk. Using the large-aperture exozodiacal spectra, we 
can still compare the shapes of the silicate features: the shape of the silicate feature from (3 Pic is 
quite different from that of comet Hale-Bopp, and it is somewhat different from that of the 
zodiacal light. While the zodiacal feature is relatively well-defined and boxy, the integrated 
exozodiacal feature is more broad and triangular. The cometary feature has a pronounced 11.3 
jjLm peak, which is not present in the zodiacal or exozodiacal spectra. Silicates producing the 11.3 
yUm feature are likely present in both exozodiacal disks; they explain the boxy shape of the red 
side of the feature. In (3 Pic, and 51 Oph, no substructure is evident in the integrated spectrum of 
the entire disk, at the limit of ISO sensitivity. The high-angular-resolution spectra (Knacke et al. 
1993; Weinberger et al. 2002) of the inner portion of the (3 Pic disk are actually similar to that of 
the zodiacal light. Any differences between the zodiacal and integrated exozodiacal may be due to 
different mineralogy, radiation environment, or mineralogy and temperatures contributing to the 
total flux. Indeed, the particle sizes and properties should be very different in the younger 
systems. While the younger systems retain their protoplanetary disks of dust and gas, the Solar 
System particles are produced by crumbling of asteroids and comets. The presence of H2O, CO, 
and CO2 disks around both 51 Oph (van den Ancker et al. 2001) and (3 Pic (Fig. 5) may lead to a 
significant gas-grain chemistry that changes particle properties. Also, while the young systems 
have very dense disks completely dominated by collisions, the Solar System is rather rarefied, and 
the dynamics of particles that produce the zodiacal light are dominated by Poynting-Robertson 



drag. 



7. Conclusions 

We have presented and contrasted the mid-infrared (5-16 fim) spectra of the zodiacal light and 

the exozodiacal light around P Pic. These two phenomena are due to planetary systems in very 
different stages of evolution: while the (3 Pic is likely in an early, particle-agglomeration stage 
leading to planet formation, the Solar System is in an advanced, large-body-colliding stage. The 
resulting particle properties, as seen in their mid-infrared spectra, are very different. For both 
systems, the mid-infrared emission is due to a combination of continuum emission and a broad 
silicate feature. The most glaring difference between /3 Pic and the Solar System is the presence 
of significant cold dust > 30 AU from /3 Pic that dominates the far-infrared emission and 
contributes significantly even in the mid-infrared. The presence of a silicate feature in both 
zodiacal and exozodiacal spectra clearly illustrates their differences. The zodiacal feature is weak, 
at only 5-7 % of the continuum, while the exozodiacal feature is brighter than the continuum. 
The amplitude of the silicate feature is an indicator of particle size. Theoretical models showed 
that the spectrum is dominated by emission from particles with radii in the range 10-100 /Ltm. 
The exozodiacal feature, in contrast, is produced by smaller particles with radii smaller than 
X/2tt = 2 /7,m. The shape of the zodiacal silicate feature, similar to the feature seen in comets and 
interplanetary dust particles, spans the range 9-12 /xm and is boxy. The exozodiacal (/3 Pic and 
51 Oph) feature is broader and more triangular. Neither the zodiacal nor the exozodiacal silicate 
features are identical to the interstellar dust feature, which peaks on the blue side and tails off on 
the red side. There is a suggestive progression in the silicate feature shape as a function of 
evolutionary age: the initial interstellar building blocks transform into the cometary spectrum 
with the P Pic spectrum as an intermediate step. 

This work is based on observations with ISO, an ESA project with instruments funded by ESA 
Member States (especially the PI countries: Prance, Germany, the Netherlands and the United 
Kingdom) with the participation of ISAS and NASA. We gratefully acknowledge help from Joris 
Blommacrt with ISOCAM calibration issues. The search for variations in the zodiacal spectrum 

was inspired by Alain Lcger. The research described in this paper was carried out at the 
California Institute of Technology under a contract with the National Aeronautics and Space 
Administration. 
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